Fractional periodicity of persistent current in coupled quantum rings 
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We study the transmission properties of a few-site Hubbard rings with up to second-nearest neigh- 
bor coupling embedded to a ring-shaped lead using exact diagonalization. The approach captures 
all the correlation effects and enables us to include interactions both in the ring and in the ring- 
shaped lead, and study on an equal footing weak and strong coupling between the ring and the 
lead as well as asymmetry. In the weakly coupled case, we find fractional periodicity at all electron 
fillings at sufficiently high Hubbard U, similar to isolated rings. For strongly coupled rings, on the 
contrary, fractional periodicity is only observed at sufficiently large negative gate voltages and high 
interaction strengths. This is explained by the formation of a bound correlated state in the ring 
that is effectively weakly coupled to the lead. 



I. INTRODUCTION 

Systems having a connected geometry exhibit inter- 
esting properties in the presence of a magnetic flux. In 
the Aharonov-Bohm effect, transmission through a ring 
pierced by a magnetic flux shows interference effects as 
electrons traversing through the two different paths ac- 
quire different phase factors. On the other hand, a mag- 
netic field induces a non-dissipativc current called the 
persistent current in isolated phase-coherent rings or ring 
ensembles due to the singlc-valucdness of the electron 
wave function— Both before and after the experimental 
observations of the persistent currents,— ~— intense the- 
oretical study has been devoted to the different factors 
affecting the periodicity of the current, such as the ef- 
fects of electron-electron interaction, disorder and finite 
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The persistent currents in more com- 



plex systems, such as multichannel rings or cylinder o 24 ' 25 
and two-dimensional quantum dot arrays? 1 - have also 
been the topic of theoretical considerations. 

The understanding of these phenomena is not only in- 
teresting for fundamental science but also from an ap- 
plicational point of view. Quantum rings could serve 
as components for future nanoelectronics. For instance, 
there have been suggestions to use ring molecules as 
quantum interference effect transistors £§r2& Aharonov- 
Bohm phenomena are also relevant for carbon nanotubes 
in which band gap oscillations as a function of external 
magnetic flux have been observed . 29 ' 30 

Metallic rings with a large number of charge carriers 
are well described by non-interacting theories and by now 
the properties of the persistent current in them, such 
as the magnitude, periodicity, and direction, are well 
understood— In these systems, the measured periodicity 
of the persistent current is a flux quantum, $ = ch/e, or 
half of it, $0/2^-^- Theoretically, the two different pe- 
riods are related to parity effects with respect to number 
of electrons in the ringi 17 ' 21 In the experiment, the half- 
flux periodicity is observed only in measurements over an 
ensemble of metallic rings.— This has been explained 
in terms of ensemble averaging of the different parities of 
the individual rings— 2 - 



In semiconducting rings with a small number of elec- 
trons, electron-electron correlations are important. With 
a strong enough electron-electron interaction, the ap- 
pearance of fractional periodicity $>o/N e i, where N e i is 
the number of electrons in the ring, has been theoreti- 
cally predicted! 16 ' 19 Fractional periodicity has also been 
shown to occur at low values of U in dilute systems when 
N e i/NU is small, associated with the disappearance of 
parity effects— 1 - However, only Keyser et aim - have been 
able to show this experimentally by studying a semicon- 
ducting ring with less than ten electrons. Very recently, 
Hernandez et aZ— - reported measurements on a 20-40- 
electron ring but they were only able to observe the $0 
and $0/2 oscillations. As the effect seems to be elusive 
in experiments, it is necessary to study which additional 
factors could complicate the experimental observation of 
fractional periodicity. 

In the experiments, either isolated rings or rings 
connected to external leads are studied. A super- 
conducting quantum interference device (SQUID)-type 
setup allows the measurement of the persistent current 
from the induced magnetic moment in isolated rings 
or ring ensembles— Also resonator-based methods 

have been used in the literature to study isolated ring 
systems— ^ These methods provide access to both the 
magnitude of the persistent current and its flux periodic- 
ity. On the other hand, the magnetoconductance of the 
ring can be measured in an Aharonov-Bohm -type setup 
with external leads— i 9 - -11 ' 33 ' 35 In these measurements, a 
bias voltage is applied, and the strength of the coupling 
between the ring and the leads can be controlled by gate 
voltages. These measurements yield only the magnetic 
flux periodicity through the measured conductance but 
provide no direct information on the magnitude of the 
persistent current. 

Theoretical study of conductance through interacting 
systems is far from trivial. A Landauer-Biittiker— -type 
formalism is available also for the interacting case, al- 
though the practical implementation is cumbersome— - 
Flux-pierced interacting rings have been widely studied 
in the literature using perturbativc approaches with re- 
spect to the coupling between the ring and the loads^^r— 
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using a variety of theoretical models to describe the ring 
part of the system: a Luttinger liquid^ spinless interact- 
ing fermions,— the t — J model,— and also the Hubbard 
models The strongly coupled case has, to the best of 
our knowledge, not yet been addressed for interacting 
systems, and in general only non-interacting leads have 
this far been considered. Moreover, we are not aware of 
any studies dealing with the conductance or transmission 
properties of Hubbard rings with second-nearest neigh- 
bor coupling. 

It has been proposed that the conductance properties 
of interacting nanostructures can be studied using the so- 
called embedding approach , 43 ! 44 in which the interacting 
system is embedded into a non-interacting lead. Instead 
of using semi-infinite leads and introducing a bias voltage 
between them, as in the widely-used Landauer-Biittikcr 
approach,— the lead is formed into a ring that is pierced 
by a magnetic flux. As a consequence, a persistent cur- 
rent is induced in the ring and its magnitude in the pres- 
ence and in the absence of the nanostructure can be re- 
lated to the transmission coefficient of the nanostructure, 
and thus to conductance. In the limit of an infinitely long 
lead, the approach yields linear response conductance. 
If electron-electron interaction is introduced to the lead 
part of the system, no formal relationship between the 
persistent current and the conductance has been found 
in the literature. However, the magnitude of the lead 
current can still be thought to probe the transmission 
properties of the nanostructure. The advantage of the 
embedding approach is that all coupling strengths can 
be treated on equal footing, and that the ground-state 
energy of the composite system is, in principle, all that 
is needed to extract the conductance properties. 

In this paper, we study a fully interacting system con- 
sisting of a flux-pierced Aharonov-Bohm ring embedded 
in an interacting ring-shaped lead. By studying a Hub- 
bard ring with second-nearest neighbor couplings, we 
show that in the weakly coupled case, the fractional pe- 
riodicity of the persistent current can be observed for all 
electron fillings. For strong coupling between the ring 
and the lead, on the contrary, fractional periodicity is 
only observable at higher electron fillings, stronger inter- 
action and larger gate voltages, where a correlated bound 
state that is weakly coupled to the lead is formed. We 
also discuss the effect of asymmetry, both in the coupling 
strengths and lead positions as well as in the interaction 
strength in the ring and lead parts of the system. Our 
results provide insight into why the fractional periodicity 
is difficult to observe experimentally. 



II. MODEL AND METHOD 

We use the Hubbard model to describe our system 
consisting of an interacting ring (hereafter ring) and an 
equally interacting ring-shaped lead (hereafter lead) con- 
nected to it, schematically shown in Fig. Q] We note 
that a related multiring or multiarm setup has recently 




FIG. 1. (Color online) The geometry of the double ring sys- 
tem. Dotted lines: alternative connections between the ring 
and the lead (the three colors/shades of gray show the differ- 
ent positions) Dashed lines: second-nearest neighbor hoppings 
in the ring. For symbols, see Section ITfl 

been studied within the non-interacting tight-binding 
formalism^ We take the hopping amplitude between the 
sites within the ring and within the lead as a constant, 
to = —t, and measure all energies in the units of t. For 
the ring part of the system, we consider rings with the 
next-nearest neighbour coupling (2NN) that is given by 
the hopping amplitude t' . 

In INN Hubbard rings, fractional §o/N e i periodicity 
appears in the U — > oo limit.— Isolated 2NN rings have 
been shown to map to continuum rings, and the inclusion 
of 2NN hoppings has been shown to lower the required 
value of the on-site electron-electron interaction U for 
the onset of fractional periodicity^ The presence of the 
2NN hopping makes thus the ring effectively quasi- ID as 
electrons are allowed to pass each other even at very high 
interaction strengths. We denote the onset value by U c . 
We use the value t' = 0.2t, following a previous study on 
isolated rings ^ for which this value was found to lead 
to a relatively low value for U c (U c = 5.3i for an eight- 
site ring at quarter-filling). The ring is pierced by an 
Aharonov-Bohm (AB) flux $ and the system properties 
as a function of this flux are studied. N r - nlg gives the 
number of sites in the ring and Ni the number of sites in 
the lead. 

In order to drive a current through the ring to probe 
the transmission properties, another AB flux that is fixed 
to the value (/>o pierces the lead ring. Throughout this 
paper, we set h = c = e = 1 and the fluxes are thus 
represented in units of $o — ch/e ~ lix. In the limit of 
an infinitely long noninteracting lead, the transmission 
through the nanostructure can be expressed using the 
current in the lead part in the presence and in the absence 
of the nanostructure, when the boundary condition in 
the lead part is given by a flux 0o = 0.25$o = tt/2, also 
called a " twisted boundary condition" i 43 i 44 i 48 We choose 
this value of 4>q also for the interacting leads, although 
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we can not analytically link the value of the current and 
the transmission coefficient. The current in the lead is 
thus driven by a hopping t e l ^° between one pair of sites 
in the lead, the rest of the hoppings between the nearest 
neighbors being to- 

The Hubbard Hamiltonian is given by 

H = i?kin + Hjj + H g , (1) 
where the kinetic energy contribution is given by 

#km = ]T(%4<V +**A C *)' ( 2 ) 

the contribution due to the on-site interaction by 

H v = UjUiin^, (3) 

i 

and the gate voltage applied upon the ring by 

ia 

Here, t%j is the spin-independent hopping amplitude be- 
tween sites i and j, Ui is the interaction strength at site 
2, and Vf a is the gate voltage applied onto the ring sites, 
modelled by an on-site energy. Other notation is stan- 
dard second quantization. We denote by U r and Ui the 
interaction strength in the ring and in the lead, respec- 
tively, if necessary. The hopping amplitude within the 
ring is given by t$e % ^ between nearest neighbours, and 
by t'e 2 ^ between second-nearest neighbours for clock- 
wise rotation.— Due to hermiticity of the Hamiltonian, 
tji = t*j naturally applies. The total flux piercing the 
ring is given by $ = N r i ng (j). The coupling between 
the lead and the ring is given by t c < 0. We assume 
there is no magnetic field present in the electron paths, 
and thus neglet the spin-orbit interaction and Zeeman 
splitting. As for a given total spin S, the states with 
different S z values are degenerate in the absence of Zee- 
man splitting, we may choose S z = (S z = 1/2) for 
even(odd)-electron systems without losing information 
on states with a higher total spin. 

We exactly diagonalize the Schrodinger equation, 
= E^, presented in the many-body basis formed 
from the many-electron states for both spin species. The 
diagonalization is performed in a block of the Hamil- 
tonian with a fixed number of electrons in the system, 
N ei = iV t + N±, and a fixed S z = iV t - N±, as the Hub- 
bard interaction does not couple sectors of the Hamil- 
tonian with different S z or N e i. To obtain the lowest- 
energy eigenstate, we apply the Lanczos diagonalization 
algorithm^, converging the lowest eigenenergy to within 
10 -12 i. The computational effort associated with the ex- 
act diagonalization method limits the number of sites in 
our system to below 30, and also the number of electrons 
in the system to around six, as the size of the Hilbert 
space grows factorially as ) ) . We mostly focus 
on a six-site ring to allow sufficient variation in the lead 



length for the consideration of finite-size effects. We also 
consider different coupling positions between the ring and 
the lead, shown by the dotted lines with different col- 
ors/shades of gray in Fig. [1] 

From the wave function of lowest-energy eigenstate, we 
calculate the observables of interest, such as the current 
between any two sites i and j, by calculating the expecta- 
tion value Ji j = (^>| Jij | ^/), where the current operator 
is given by (Ti = c = e = 1)££ 

Jij = * ^](^ i J C L c JO' — ^ij' c ](T C io-)- (5) 
a 

We define J r i ng and Ji ea d such that the sites i and j in 
Eq. ([5]) both belong to the ring or lead part of the sys- 
tem, respectively. Alternatively, we could calculate the 
persistent current in the ring from the variation of the 
ground-state energy E with respect to a change in the 
fluxi£ 



We prefer, however, the current operator approach as it 
allows us to easily calculate the current in the different 
branches of the system by choosing the sites i and j ac- 
cordingly. Similarly, we calculate the occupation of a 
given site (hi = ^Z a h la = J2a c L c ^) or tlic total s P m 
of our system [(S 2 ) = S to t(Stot + 1)] by calculating the 
expectation value of the corresponding operator. 

III. RESULTS 
A. Weakly coupled rings 

To begin with, we consider weak coupling between the 
ring and the lead. We expect the behavior of the frac- 
tional periodicity to be close to that of isolated rings, and 
also to observe single-electron tunneling. For the time 
being, we restrict ourselves to the symmetric ring-lead 
coupling (Fig. [I] 1-4 connection). 

We use a 14-site lead and leave the discussion on finite- 
size effects to Section UTTCl We occupy the system with 
an odd number of electrons, — 2 and = 3. The 
ring and the lead are weakly coupled, t c = —0.2t. We 
calculate the current in the ring and lead parts of the 
system, as well as the total spin state of the composite 
system and its contribution from the ring, at V 9 values 
chosen as to keep the ring occupation n r ; ng near an in- 
teger value, as shown in the top left panels in Fig. [2] at 
U = 2t, 5t, lit, and 21t. The values of U have been cho- 
sen to best illustrate the fractional peridicity §o/N e i for 
each of the ring fillings N e i = 2,3,4 and 5. To restrict the 
computational effort, however, only integer values for U 
were considered. In the right panels of Fig. [^a)-(d), the 
current in the ring, J r i ng , is given by the solid lines and 
the one in the lead, Ji oa d, by the dotted lines. The cur- 
rent have been shifted and normalized to better illustrate 
the shape and periodicity of the ring current. 
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FIG. 2. (Color online) The properties of a six-site 2NN ring weakly connected (t c = — 0.2t) to a 14-site lead as a function of 
the total flux $ piercing the ring, with the ring occupation 72 r j n g at integer values (top left panels). Bottom left panels show 
the total spin (thin line) and its projection onto the ring part of the system (thick line) and the right panels the ring current 
(solid lines) and the lead current (dashed lines). The dotted lines show the zero level, and the curves have been normalized to 
facilitate comparison. For lead currents, the scaling factors rounded to two significant digits are indicated. The values of U 
are chosen as to show fractional periodicity at each integer n r j ng , and the gate voltage V 3 to keep the occupation constant as 
a function of the flux $. The five curves illustrate the different ring occupations (n r i ng « 1,2,3,4,5: black, blue, red, green, 
yellow in color, black to light gray in grayscale), (a) U — 2t (b) U = 5t (c) U = lit (d) U = 21t. 
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FIG. 3. The two lowest energy levels as a function of the 
flux $ for a six-site ring with Ni = 14 as given by the Lanc- 
zos diagonalization in which the lowest-energy state has been 
converged. The ring occupation n r i ng is approximately three 
electrons, (a) Weakly coupled system at fractional periodic- 
ity (U = 5t, t c = -0.2t, V 9 = -0.9t). (b) Weakly coupled 
system for U > U c . (U = St, t c = -0.2t,V 9 = -l.lt) (c) 
Strongly coupled system near fractional periodicity (17 = 8t, 
t c = -t, V 9 = -l.lt). Other parameters tt = 0.2t, N el = 5. 



In Fig. [2] (right panels, solid lines), the kinks in the 
ring current appear due to crossings between the energy 
parabola that determine the ground state of the sys- 
tem. Fig. [3] illustrates this, showing the energy of the 
two lowest-energy states as a function of i> as obtained 
from the Lanczos diagonalization in which the lowest- 
energy state has been converged. The energy parabo- 
las are easily identified. At the crossings, the slope of 
the ground-state energy changes sign, leading also to an 
abrupt sign change in the persistent current [see Eq. flSJ)]. 
For n r i ng = 1, the periodicity of the ring current is always 
$o- For higher occupations, in the absence of interactions 
a period of $o ($o/2) would be expected for even (odd) 
occupation numbers in the ringi 16 ' 17 Interaction, how- 
ever, shifts the energy parabola, leading to additional 
crossings and a pseudo-$o/2 period for even electrons 
below ET ^ 16 ' 17 as parabola corresponding to states with a 
higher total spin are less shifted. This is also reflected 
in the total spin projected onto the ring, shown in left 
panels in Fig. [2] by thicker lines. When U is increased, 
the regions where the ring spin is high or even maximal 
(S = riring/2) are widened. 



At the onset of the fractional periodicity (U = U c ), 
the period of the ring current is &o/n T i ng and segments 
between the kinks are of equal length in $.This is also il- 
lustrated in Fig. EHa) that shows the ground-state energy 
as a function of <£> at U = 5t corresponding to n r j ng = 3 
with fractional periodicity. Notably, even though two of 
the ring sites are perturbed through the presence of the 
leads, the fractional periodicity is not affected. Above U c , 
the parabola continue to shift in energy and the kinks in 
the ring current are no longer of equal length, shown in 
Fig. [3Jb) at U = 8t. It is readily seen that the parabola 
around $ = $o/2 + n$o have shifted relative to the ones 
at $ = ri$o- The periodicity, however, remains pseudo- 

$o/«rmg- 

The different periodicities of the persistent current in 
the ring are also reflected in the lead current, shown in 
dashed lines in the right panels of Fig. [2] Apart from the 
Turing — — 5 case, in which the magnitude of the ring current 
is so small that it is significantly affected by numerical 
instabilities, the kinks in the ring current are associated 
with jumps or peaks in the lead current. The overall 
periodicity of the lead current is twice that of the ring 
current. This is understood by considering the double- 
ring geometry (Fig. [I}. In the present configuration (1-4 
connection, see Fig. [I}, the ring is pierced by the flux 6(f>. 
The lead and some of the ring sites can be thought to 
form another ring that is pierced by the flux — 30 + 4>q. 
Thus, the periodicity in the lead current is twice that of 
the ring current. The effect of the position of the coupling 
sites will be discussed in more detail in Section MI Dl 

In Fig. [2J the lower left panels show the total spin of 
the ring-lead system (thin line) and its projection onto 
the ring part (thick lines). The total spin is restricted 
to half-integer values by the electron number (n-j- = 2, 
ni = 3), and at low interaction values it is 1/2 for all 
"ring- At higher U, some regions of S = 3/2 appear. 
The total spin projected onto the ring, on the contrary, 
reflects the changes in the ground-state energy parabola. 
With an increasing interaction strength, the regions with 
a higher total spin value widen in $, and for n r i ng = 4 
even regions of maximal spin 5 = 2 appear at and above 
U c . This effect is similar than in isolated rings^i in which 
a maximally spin polarized state (S = N e i/2) was found 
above U c at a fixed value of $ = 0.5$o for a quarter-filled 
eight-site ring. Due to the coupling between the lead and 
the ring, the values of the total spin projected onto the 
ring are, however, not strictly integer or half-integer. 

The onset value of the fractional periodicity scales ap- 
proximately as U c oc N^i- Thus, in order to experi- 
mentally observe fractional periodicity, a large number 
of electrons in the ring requires very strong interaction 
that might not be experimentally realistic. In the weakly 
coupled case, the fractional periodicity is very similar to 
the isolated case and we thus turn our attention to the 
case with a strong coupling between the lead and the 
ring. 
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B. Strongly coupled systems 

Similar to Fig. [5] for weak coupling, Fig. @] shows the 
ring and lead currents, along with the total spin and the 
ring occupation in the strongly coupled systems. Again, 
the values of U are chosen as to best show the fractional 
periodicity at each n r j ng if possible. The top left panels 
in Fig. 0] show the ring occupation as a function of $. 
Unlike for weak coupling, the charge state of the ring can 
not be fixed to integer values for all electron numbers and 
interaction strengths. For n r i ng = 2, the value U = 2t was 
chosen to facilitate comparison with the weakly coupled 
case (Fig. [5]), and for n ring = 3, U = 8t to show the 
closest the system gets to fractional periodicity. 

The ring and lead currents (Fig. 01 right panels) for low 
occupations appear irregular and in the n r i n g = 1 case 
even 2<& periodic. This is, again, related to the doubled 
period in the lead part of the system. In the strongly 
coupled case, the magnitudes of the lead and ring cur- 
rent are comparable and due to the current conservation 
at the coupling points, the lead current alters the ring 
current. In contrast to the weakly coupled case, we only 
observe fractional periodicity for higher ring occupations, 
"ring > 4. For electron fillings n r in g < 3, the saw-tooth 
shape of the persistent current is completely smoothened 
at low interaction strengths, and for n r j ng = 3 a persistent 
current resembling fractional periodicity with a partly 
smoothened shape appears. This is related to entangle- 
ment with the lead states and the resulting fluctuation in 
the occupation of the ring as a function of $ i 14 i 25 On the 
other hand, from the viewpoint of the ring, the strong 
coupling to the lead can be interpreted as a form of dis- 
order at the coupling sites. Disorder has also been shown 
to lead to smoothening of the kinks^ In the n r i ng = 3 
case, however, we observe a ring current characteristic to 
U > U c in the weakly coupled case at strong interaction, 
even though the actual fractional periodicity with sharp 
saw-tooth features is missing. 

At higher ring occupations, n r j ng = 4 or 5, we see frac- 
tional periodicity at sufficient U, and (pseudo-)$o/2 pe- 
riodicity below it. Thus, it appears that at sufficiently 
low V 9 values corresponding to high n r ing) the ring be- 
comes effectively weakly coupled to the lead. Some of the 
kinks are slightly smoothened due to avoided crossings of 
the energy parabola, fluctuations in the ring occupation, 
and large magnitude of the lead current. Fig.[3](c) shows 
the ground-state energy as a function of $ at U = 8t 
for n r i ng w 3. As opposed to the weakly coupled case 
of Fig. [3(b), some of the sharp crossings between the 
parabola disappear in the strongly coupled case. For oc- 
cupations and interaction strengths with fractional peri- 
odicity, the behavior of the lead current, as well as the 
behavior above U c , is similar to that at weak coupling. 
For values of V 9 for which the ring occupation strongly 
fluctuates as a function of $, also other periodicities are 
possible. For instance, in a six-site ring strongly cou- 
pled to a 14-site lead with five electrons, at V 9 ~ — 1.5i 
and U = 2t the occupation fluctuates between three and 



five electrons, and the periodicity of the ring current is 
pseudo-$ /3 (not shown). 

In general, in the strongly coupled systems the peaks 
and kinks in the lead current do not reflect the kinks in 
the ring current as accurately as for weak coupling. At 
fractional periodicity, the behavior of the ring and lead 
currents arc similar irrespective of the coupling strength, 
namely, with abrupt magnitude and sign changes that 
coincide with the jumps in the ring current. At low U 
and lower ring occupations we mainly see a single peak 
in the 2$o _ P er iodic lead current, as if the ring were just 
a disordered site in the lead ring. For n r i ng = 5, the lead 
current shown in Fig. U is essentially numerical noise due 
to minuscule magnitude of the current. Thus, in the 
strongly coupled case, the observation of fractional peri- 
odicity from conductance or transmission is much more 
difficult than at weak coupling. 

The evolution of the total spin with increasing interac- 
tion is similar to the weakly coupled case, although the 
spin projection onto the ring is smoothened out at low 
"ring due to the strong coupling. At fractional periodic- 
ity, however, the crossings between different- S parabola 
are clearly seen as abrupt changes in the spin state in the 
upper right panels of Fig. 2J Again, with increasing U, 
states with a higher total spin become the ground state 
for larger regions in $. 

In Figs. [2] and [4l the current curves have been normal- 
ized as to facilitate qualitative comparison. It is, how- 
ever, instructive to also consider the magnitudes of the 
current. For the lead currents, the scaling factors relative 
to the n r i ng = 1 current are given in the figures, rounded 
to two significant digits. In general, the magnitude of 
the ring current is only slightly smaller in the strongly 
coupled case than in the weakly coupled case, the dif- 
ference increasing with increasing interaction strength. 
The lead current, on the other hand, is an order of mag- 
nitude larger in the strongly coupled case than in the 
weakly coupled case. In particular, for the strongly cou- 
pled n r ; ng = 1 case, the lead and ring current are of 
the same order of magnitude. The strong coupling dis- 
turbs the ring part, decreasing the ring current as the 
hybridization with the lead states is enhanced. The in- 
creased magnitude of the lead current in the strongly 
coupled case can be understood as increased transmis- 
sion through the ring. 



C. Finite-size effects 

As we are dealing with relatively small systems, it is 
necessary to pay attention to finite size effects, namely, 
what is the effect of the length of the lead. We are not in- 
terested in the diminishing magnitude of the lead current 
with increasing lead length^ but instead, the behavior of 
the ring persistent current and the gate voltages at which 
the charge state of the ring changes. Fig. [5] shows the 
properties a six-site ring with 2NN hopping (t 1 = 0.2t) 
with the length of the lead increasing from 10 sites to 20 
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FIG. 4. (Color online) The properties of a six-site 2NN ring strongly connected (t c = —t) to a 14-site lead as a function of the 
total flux $ piercing the ring, with the ring occupation n r i ng as close as possible to integer values (top left panels). Bottom left 
panels show the total spin (thin line) and its projection onto the ring part of the system (thick line), and the right panels the 
ring current (solid lines) and the lead current (dashed lines). The dotted lines show the zero level, and the curves have been 
normalized to facilitate comparison. For the lead currents, the scaling factor rounded to two significant digits is indicated. The 
values of U chosen as to show fractional periodicity at each n r i ng , or as close as the system gets to it, and the gate voltage V 9 
to keep the occupation constant as well as possible. Due to the strong coupling, large fluctuations in the electron number are 
present especially at low values of interaction. The five curves illustrate different ring occupations (n r i ng ~ 1,2,3,4,5: black, 
blue, red, green, yellow in color, black to light gray in grayscale), (a) U = 2t (b) U = 8t (c) U = 16t (d) U = 22t. 
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FIG. 6. (Color online) The effect of an increasing lead length 
Ni on the Coulomb staircase in a six-site ring strongly coupled 
(t c — —t) with the lead. In (a), N e i = 5, the lead lengths 10, 

13, 16, and 19 sites, and in (b) N e i = 6, the lead lengths 11, 

14, and 17 sites, the lead length increasing with darkening 
shade of gray/red. In (c), the Coulomb staircase for different 
electron numbers is compared (N e i —5: Ni — 19, N e i = 6: 
Ni = 17, and N el =7: JVj = 14). In all calculations, = 
0.01$o and U = 2t. 



FIG. 5. The effect of lead length at weak coupling (t c — 
— 0.2t) with a fixed flux <f>i = O.Ol'I'o on (a) the occupation of 
the ring (b) the lead current as a function of the gate volt- 
age V 9 . (c) The persistent current in the ring (d) persistent 
current in the lead as a function of the flux piercing the ring, 
at V g chosen such that n r i ng = 2. The number of sites in 
the lead, Ni, increases from 10 to 20 in steps of two with the 
darkening shade of gray. Other parameters N e i = 6, U = 2t 
and t' = 0.2t. 



sites in steps of two sites. The ring is weakly coupled, 
t c = —0.1 i, and the number of spin up- and spin down- 
electrons are equal (N e i = 6). In (a) and (b), fa is fixed 
to 0.01$o and in (c) and (d), V 9 in the ring is chosen to 
give n ri n g ~ 2. 

In Fig. [5ja), a clear Coulomb staircase characteris- 
tic for weakly coupled systems with well-defined charge 
states and transition regions associated with an electron 
tunneling from the lead to the ring between them is seen. 
As the lead length is increased, the steps shift toward 
more negative gate voltages, and the values of V 9 cor- 
responding to the steps converge starting from the steps 



with highest n r ; n g- As the total number of electrons in 
the ring-lead system is fixed, a removal of an electron 
in the ring corresponds to an addition of an electron in 
the lead. For high ring occupations, the filling factor 
in the lead is low and thus the exact number of sites is 
not as important as for higher fillings. The lead can be 
thought to act as a sort of a small reservoir but due to the 
system configuration, neither the electron number of the 
ring nor its chemical potential remain constant when the 
gate voltage or flux piercing the ring are changed. Com- 
paring the lead current as a function of the gate voltage, 
shown in Fig. [5jb) , we see that changes in the ring occu- 
pation are accompanied by peaks in the lead current due 
to charge transfer between the two subsystems. 

Figs. [5jc) and [5jd) show the persistent current in the 
ring and in the lead as a function of $, respectively, at a 
gate voltage V 9 fixed such that there are approximately 
two electrons in the ring. The small differences between 
the curves in Fig. (5jc) are mainly due to small fluctua- 
tions in the ring occupation. We see that the lead length 
has little effect on the current within the ring, whereas 
with an increase in the lead length leads, the magnitude 
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FIG. 8. The effect of asymmetry in the strength of the lead- 
ring coupling. The product of the couplings is fixed to one, 
tcit C 2 = 1, but the ratio between them is varied. Parameters: 
Six-site ring, Ni = 14, U = 5t, N et = 5, fa = 0.01$o, and 
t' = 0.2t. 



of the lead current is somewhat decreased [Fig. EJd)]. 
The persistent current in a perfect, non-interacting ring 
without a nanostructure decays as l/Mead, and with a 
weakly coupled embedded Anderson impurity-type quan- 
tum dot as l/A^old-" Again, qualitatively the curves are 
very similar, the changes in the current occurring at same 
flux values with sign changes in the ring current. 



In the strongly coupled case, the Coulomb staircase 
is partly smoothcncd out, as illustrated in Fig. |6] In 
principle, the strong coupling between the lead and the 
ring should lead to broadening of the transmission reso- 
nances leading to the smoothening of the Coulomb steps. 
Fig. H2a) and (b) show the ring occupation as a function 
of the gate voltage V 9 for a six-site 2NN ring at strong 
coupling (t c = — t) and U = It for N e i = 5 and N e i = 6, 
respectively. In both subfigures, the darkening shade of 
the curves indicates the increasing number of sites, Ni, 
in the lead. We indeed observe a smoothening of the 
steps, to the extent that in the N e i = 5 case only one 
step is observable, compared to all five seen at weak cou- 
pling. Fig. EJc) compares the staircase for N e i = 5, 6, 
and 7 with Ni = 19, 17, and 14, respectively. The low- 
occupation steps remain smoothcncd out but increasing 
the total number of electrons in the system leads to the 
presence of additional high-occupation steps. 

The disappearance of the low-occupation steps can be 
understood if the competition between kinetic and inter- 
action energy is considered. At large negative gate volt- 
ages, the localization of electrons in the ring is favored 
as the kinetic effects due to the strong coupling are small 
compared to the effect of the gate voltage. As the ring is 
close to half-filling, the density of doubly occupied sites 
is low and due to the presence of the 2NN hoppings, the 
electrons are allowed to move within the ring without 
having to doubly occupy any single site. When the gate 
voltage is increased, the electrons are increasingly delo- 
calizcd both in the ring and the lead. Charging effects 
associated with the addition of a single electrons can no 
longer be observed. Thus, at large negative V g the sys- 
tem is effectively weakly coupled to the gate voltage, and 
shows behavior similar to weak coupling. At higher V g , 
on the contrary, neither well-defined charge states in the 
ring, nor fractional periodicity in the ring current are 
observed. 

Studying the effect of the interaction strength U on the 
Coulomb staircase provides additional insight into this 
effective decoupling. Fig. [7]Ja) shows the effect of the 
interaction strength on the Coulomb steps. When the 
interaction is increased, shown in the different shades of 
gray, the number of steps increases. With an increas- 
ing electron-electron repulsion strength, the kinetic en- 
ergy associated with the strong coupling becomes weaker 
relative to the interaction energy, and bound correlated 
states are formed in the ring part at sufficiently negative 
V 9 also for lower n r i ng . The bound correlated states also 
support a ring current with fractional periodicity In- 
deed, we see in Fig. [Tfa) that at U > 15i, a Coulomb 
step with n r i n g ~ 3 has appeared, and the corresponding 
n ring = 3 curves in Fig. U show characteristics increas- 
ingly similar to the weakly coupled case above U c with 
increasing U. 

The interaction strength in the lead is actually less im- 
portant. Figs.JTJb) and [7|c) show the effect of increasing 
the interaction in the lead part of the system, while the 
interaction in the ring part is U r = 5t and U r = 20t, 
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FIG. 9. (Color online) The ring and lead currents [(a)-(c) and (d)-(f), respectively], as well as spin imbalance — n^ in the ring 
[(g)- (i)], in the (<£, V 9 ) space for the three different coupling sites for weak coupling between the lead and the ring. (a),(d),(g): 
1-2 coupling (b),(e),(h): 1-3 coupling (c),(f),(i): 1-4 coupling. In (a)-(c), the contours of half- integer ring occupation n r i ng are 
additionally drawn, separating regions of approximately integer values of n r i ng . Other parameters: Six-site ring, 14-site lead, 
t c = — 0.2t, t = 0.2i, N e i = 5, U = 2t. The color scales for each row are fixed to allow comparison in the magnitudes between 
different coupling positions. 



respectively, for all curves. The highest- V 9 steps slightly 
shift when the interaction is increased but the overall 
shape of the curve is very similar regardless of the lead 
interaction strength. At weaker interaction in the ring 
(U r = bt in Fig. [TJb)) the presence of interaction in the 
lead has little effect, whereas of stronger ring interaction 
(U r = 20t, Fig.[7][c)), the system is practically converged 
to the Ui = 20t case at Ui = bt. 



D. Asymmetry effects 

This far, we have only considered structures, in which 
the coupling between the lead and the ring is symmet- 
ric, both in strength and in the coupling position. In 
experiments^ however, the coupling strength between 
the ring and the lead is typically controlled by gates and 
it is possible that the coupling strength is not equal for 
both contacts. Our approach allows us to study explic- 
itly the effect of both types of asymmetry. Fig. [5] com- 
pares the Coulomb steps for a six-site ring 1-4 -coupled 



to a 14-site lead but with unequal coupling strengths. 
The product of the couplings, t c \t c i = 1 is fixed but the 
asymmetry, given by the fraction t c \/t C 2, is varied. When 
the difference between t c \ and t C 2 increases, the systems 
becomes effectively more weakly coupled, seen in the ap- 
pearance of the Coulomb steps. Thus, a single weak link 
is sufficient for weakly coupled behavior. 

The position of the lead-ring coupling can also be 
changed, the arms of the ring thus being of un- 
equal length. Unequal-lenght arms appear also in the 
experiments^ Fig. Q] illustrates the three different cou- 
pling positions available for the six-site ring. The effect 
of the coupling positions on the ring and lead currents 
is illustrated in Fig. [9] in which the magnitude and sign 
of the ring current [(a)-(c)] and the lead current [(d)- 
(f)] is shown in the (V 9 ,$) space for a weakly coupled 
(t c = — 0.2t) six-site ring with a 14-site lead (U = 2t, 
t' = 0.2£, n-f = 3, ni = 2), as well as the spin imbalance 
in the ring, n-\ — Hj, [(g)- (i)]. The color scale is fixed so 
that in each row, the magnitudes for different coupling 
positions can be compared. The contours overlaid with 
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FIG. 10. (Color online) The ring and lead currents [(a)-(c) and (d)-(f), respectively], as well as spin imbalance rif — nj. in 
the ring [(g)-(i)], in the ($, V 9 ) space for the three different coupling sites for strong coupling between the lead and the ring. 
(a),(d),(g): 1-2 coupling (b),(e),(h): 1-3 coupling (c),(f),(i): 1-4 coupling. In (a)-(c), the contours of integer and half-integer 
ring occupation n r i ng are additionally drawn. Other parameters: Six-site ring, 14-site lead, t c = —t, t' = 0.2t, N e i = 5, U = 2t. 
The color scales for each row are fixed to allow comparison in the magnitudes between different coupling positions. 



the ring current maps in Fig. IHta)-(c) show the occupa- 
tion of the ring, n, r i ng . 

It is easily understood that for weak coupling, the 
choice of the position has little effect on the behavior 
of the ring current and fractional periodicity, as seen in 
Fig. HJa)-(c). The effect of the position is, however, seen 
in the current circulating in the lead part of the sys- 
tem. More precisely, an alternation in the magnitude of 
the lead current depending on the coupling position is 
clearly seen, as well as changes in the broadening of the 
lead current peaks. The high-current regions approxi- 
mately parallel to the $ axis are due to the electrons 
entering the ring, seen in the upper row as contours of 
half-integer n r i ng . The ridges parallel to the V 9 axis, on 
the other hand, are due to charge transfer changing the 
spin polarization n,| — nj_ within the ring at a constant 
ring occupation. The density of one spin species increases 
and the other decreases abruptly, as seen in Fig. IUJg)-(i). 

In fact, the magnitude alternation in the lead current 
[Fig. injd)-(f)] resembles the magnitude of persistent cur- 
rent in a ring, superimposed on the resonance peaks, only 
with a two, three or six times longer period. There is 



also a phase shift when the number of electrons in the 
lead changes, similar to the odd-even alternation in iso- 
lated, clean rings. Like already mentioned, the effect of 
increased period in the lead current is seen also in the 
symmetrically coupled case. In Figs. [2] and |H the peri- 
odicity of the lead current is twice that of the ring cur- 
rent, explained by a second ring formed by the lead and 
some ring sites that is pierced by half of the flux $ in 
Section IIII Al Also for the other coupling positions, we 
get the period of the lead current by comparing the flux 
piercing the ring, and a ring formed by the lead and some 
sites in the ring. Thus, the 1-2 coupling results in a six 
times longer period in the lead current compared with the 
period of the ring, and the 1-3 coupling in a three times 
longer period. We also note that the maximal magnitude 
of the lead current in the symmetrically coupled case is 
approximately half of that for the other cases, shown in 
Fig. EIc). This is most likely due to the appearance of 
destructive interferences when the ring arms are of equal 
length. 

Fig. [TU] shows similar figures for the strongly coupled 
case. As there are no bound correlated states in the V 9 
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FIG. 11. (Color online) The occupation in the ring, Wring as 
a function of the gate voltage V 9 for ring sizes iVring = 5 — 8 
for N e i = 5, 6, and 7 (from lighter to darker curves). For 
the five-site ring, also N e i = 8 included. Other calculational 
parameters: U = 2t, t' = 0.2t, t c = -t, and iV; = 14. 



range considered and the electrons are delocalized both 
in the ring and the lead, the coupling position affects also 
the ring current. For larger negative V 9 , regions with a 
well-defined constant n r i ng similar to the weakly coupled 
case would be seen also for strong coupling. The effect 
of the coupling position on the magnitude of the lead 
current is much larger than in the weakly coupled case, 
the maximal value of the lead current being six times 
larger at 1-2 -coupling compared to 1-4 -coupling. The 
diminishing magnitude of the lead current also makes the 
charge imbalance smaller in the 1-4 case, shown in the 
contours of the lower row in Fig. 1101 The periodicities 
of the lead current, however, are the same as at weak 
coupling. 

Up to this point, we have only considered a six-site 
ring. We conclude by noting that our results on the 



six-site rings are reprcscntable also for other ring sizes. 
Fig. QT] compares the Coulomb steps at strong coupling 
for five- to eight-site rings with a 14-site lead. In all cases, 
the lowest-occupation steps are smoothened out, and for 
higher n r j ng the bound states that are effectively weakly 
coupled to the lead appear. 



IV. CONCLUSIONS 

We have studied a Hubbard ring with second-nearest 
neighbor hoppings, connected to a ring-shaped lead using 
an embedding approach and solving the system Hamil- 
tonian using exact diagonalization. In the case of a weak 
ring-lead coupling, the behavior of the persistent current 
is similar to that in isolated rings and the presence of the 
lead is only a weak perturbation. In addition, the inter- 
action strength required to observe fractional periodicity 
is increased with a higher electron occupation in the ring. 

In the strongly coupled case, on the contrary, fractional 
periodicity can not be observed for low electron occupa- 
tions in the ring. For higher occupations corresponding 
to large negative gate voltages, fractional periodicity is 
observed due to the formation of bound, strongly cor- 
related states. High interaction strength, however, are 
required. Our results provide additional insight to why 
the fractional periodicity is difficult to observe experi- 
mentally. 
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